Quantum criticality in spin chains with non-Ohmic dissipation 
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We investigate the critical behavior of a spin chain coupled to bosonic baths characterized by a 
spectral density proportional to ui s , with s > 1. Varying s changes the effective dimension d c ff = d+z 
of the system, where z is the dynamical critical exponent and the number of spatial dimensions d is 
set to one. We consider two extreme cases of clock models, namely Ising-like and [/(l)-symmetric 
ones, and find the critical exponents using Monte Carlo methods. The dynamical critical exponent 
and the anomalous scaling dimension r\ are independent of the order parameter symmetry for all 
values of s. The dynamical critical exponent varies continuously from z fa 2 for s = 1 to z = 1 for 
s = 2, and the anomalous scaling dimension evolves correspondingly from n > to rj = 1/4. The 
latter exponent values are readily understood from the effective dimensionality of the system being 
d e $ « 3 for s = 1, while for s = 2 the anomalous dimension takes the well-known exact value for the 
two-dimensional Ising and XY models, since then d c ff = 2. However, a noteworthy feature is, that 
z approaches unity and rj approaches 1/4 for values of s < 2, while naive scaling would predict the 
dissipation to become irrelevant for s = 2. Instead, we find that z = l,r] = 1/4 for s m 1.75 for both 
Ising-like and U(l) order parameter symmetry. These results lead us to conjecture that for all site- 
dissipative Z q chains, these two exponents are related by the scaling relation z = max{(2 — rj) /s, 1}. 
We also connect our results to quantum criticality in nondissipative spin chains with long-range 
spatial interactions. 
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I. INTRODUCTION 

The spin-boson modePEl (SBM) represents one of the 
most well-established frameworks for describing the effect 
of dissipation on a quantum system. In its simplest incar- 
nation, it describes a two-level system coupled to an in- 
finite number of harmonic oscillators with low-frequency 
spectral density J(oj) oc uj s , with Ohmic damping (s = 1) 
being the most commonly studied case. Generalizations 
of this model include extensions to finite (spatial) di- 
mensions d > Oj^and models where the Z 2 (Ising) s pin 
symmetry has been replaced by a higher symmetryPEl 
Extended versions of such systems may also find applica- 
tions in the study of quantum critical points in quantum 
magnets and strongly correlated systems,®^ and hence 
they are of considerable interest in contemporary con- 
densed matter physics. 

Another generalization is to consider non-Ohmic spec- 
tral densities (s ^ 1), which may be relevant in the 
description of several different physical phenomenaP^HH] 
From a more fundamental physics point of view, the sub- 
Ohmic (s < 1) SBM and related models have been stud- 
ied intensively in recent years^^Sl following claims that 
the so-called quantum-to-classical mapping may be vio- 
lated in even the simplest variant of SBM for s < 1/2 
Its super-Ohmic counterpart s > 1 has, on the other 
hand, received relatively little attention. This may be 
due to the fact that the (0+l)-dimensional [(0+l)D] SBM 
exhibits a (quantum) phase transition only for values of 
s < 1. For d > 1, however, the possibility of a phase 
transition arises for all s. 

The SBM is generally described, via the quantum-to- 
classical mapping, by a classical (d+l)D spin model with 



long-range interactions that decay as l/r 1+s in imaginary 
time t. Long-range interactions are interesting, as they 
allow one to increase the effective dimensionality con- 
tinuously by tuning s to lower values. In classical spin 
glasses, for instance, low-dimensional models with long- 
range interactions have been studied to infer properties of 
higher-dimensional realizations of the same systems with 
purely short-range interactionsJ^^HUl j n quantum models, 
the effective dimensionality is expressed by d c ff = d + z, 
with z being the dynamical critical exponent defined from 
the divergence of the correlation time £ r ~ £ z , where £ is 
the spatial correlation length. At a second-order phase 
transition, we have in standard notation £ ~ \K — if c | _l/ 
as the coupling parameter K approaches its critical value 
K c . The presence of dissipation in general causes z to de- 
viate from the value z = 1, with a naive scaling estimate 
zq = Although this result is exact in mean-field 

theory (d c g > 4), deviations may appear when decreas- 
ing d e ff. For the Ohmic case, it is knowrP^that z obeys 
the scaling law z = zq — rj, where rj in general denotes 
the anomalous scaling dimension at the transition point 
to a disordered state. Below d e s = 4 one has r\ > 0, and 
previous worlPl on s = 1 for d = 1 found 77 0.015 and 
z w 1.985. 

One issue we address in this paper is how the exponents 
z and 77 evolve as one varies the dissipation parameter 
s > 1 . For the Ohmic case considered previously, the de- 
viations from naive scaling (i.e., from z = zq) are barely 
significant due to the small value of rj when d c g ps 3. This 
deviation should become more noticeable as the effective 
dimensionality decreases, although one cannot expect the 
relation z = zq — rj, valid for s — 1, to hold also for larger 
s. In the limit d+z — > 2, the anomalous dimension 
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might be expected to approach the relatively large value 
r\ = 1/4, which it takes for both the two-dimensional (2D) 
Ising and 2D XY model. A related issue is the value of 
s beyond which the dissipation term is irrelevant in the 
renormalization group sense, giving z — 1. Naive scal- 
ing indicates that z = 1 for s > 2, but as z is likely to 
decrease faster than z = 2/s as s increases, dissipation 
might turn irrelevant for a value of s smaller than 2. 

Another issue which we address is how the critical ex- 
ponents, in particular z and ?y, depend on the symmetry 
of the order parameter. In the limit s — 1, there is no sig- 
nificant difference between the values of 77 (and ther eby 
z) for discrete and continuous order parameter fields 
but it is not inconceivable that such a difference becomes 
noticeable for lower effective dimensions, that is, as s is 
increased. 

In order to answer these questions and to study a class 
of dissipative models for which relatively little is known 
precisely, we present results from Monte Carlo simula- 
tions on both XY and Ising-like spin chains with non- 
Ohmic dissipation. In both cases, we consider super- 
Ohmic dissipation, which for the XY case allows us to 
interpolate between the universality class describing the 
three-dimensional (3D) XY model and the very differ- 
ent Berezinskii-Kosterlitz-Thouless (BKT) criticality of 
the 2D XY model. The models are presented in the 
next section, where we also describe the finite-size scal- 
ing procedure used to find the critical exponents. The 
dependence of these exponents on s are presented and 
discussed in Se c. |III| before we give a summary of our 
findings in Sec. |IV| 



II. MODEL AND FINITE-SIZE SCALING 
METHODS 

The starting point of the models we consider may be 
taken as a general (1+1)D 4 -type quantum field theory 
of an O(N) order parameter field (f>. Including dissipa- 
tion, the Fourier transform of its inverse bare propaga- 



tor is of the form q 



\u>\ s , where the damping 



term oc \ui\ s arises from the coupling of the field to baths 
of harmonic oscillators^ with a low-frequency power-law 
spectral density characterized by the exponent s. 

Parameterizing the order parameter field of such a N = 
2 quantum rotor model by an angle variable 9, we may 
formulate the discretized action as 
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on a quadratic L x L T lattice. Above, K is the spa- 
tial coupling constant to be varied, whereas the quantum 



coupling constant K T and the dissipation strength a are 
taken as fixed values during the simulations. 

In order to study both continuous and discrete symme- 
try of the order parameter field, we consider two possible 
domains of the angle variables: f7(l) symmetry is equiv- 
alent with 9 G [0,27r), and for a discrete symmetry (Z4), 
we choose to enforce the restriction 9 e {0, |-,7r, ^f}. 
We refer to the former as the XY model and to the 
latter as the model. Such a Z 4 model will be in 
the same universality class as a corresponding Z2 (Ising) 
model, which is why we refer to this model as Ising-like. 
This equivalence is easily shown using the substitution 
cos 9 X . T = (a x>T + Mz,r)/2 and sin 6^ = (<t X)T - £t x , T )/2 
to rewrite the action as that of two Ising models in 
terms of decoupled Ising spins a and /i. The order pa- 
rameter for both symmetry variants is defined as m — 
(LL T )~ l ^ exp (i9 x , T ) in the standard manner. 

When determining critical exponents of a quantum sys- 
tem using finite-size scaling (FSS), the system dimensions 
L, L T used have to be chosen such that they respect the 
system anisotropy reflected by the dynamical critical ex- 
ponent, L T oc L z . This is a problem when we do not 
know dynamical critical exponents a priori, and one usu- 
ally has to first determine z by simulating several values 
of L T for each L, before running new simulations with 
L/L* fixed. We circumvent this problem by using the 
same data to determine z and to evaluate the FSS ob- 
servables, by interpolating data for multiple L T values 
to L T = L* T (L). Here, L* is a characteristic temporal 
system size found for each spatial system size L, as ex- 
plained below, and it is assumed that L* oc L z . This has 
the advantage that one (along with z) can find all other 
critical exponents simultaneously, utilizing all (or most 
of) the data generated. Furthermore, we are also able 
to appropriately take the uncertainty in z into account 
when finding the uncertainty in the other exponents, by 
repeating the entire procedure for a number of jackknifc 
bins based on the original data. 

The procedure to find z is explained in more detail in, 
for example, Ref. 131 and is based on the Binder ratio 
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Q(\K-K C \I}I\L T IL Z ), 



(2) 



where brackets (...) indicate ensemble averages and Q 
is a universal scaling function. The characteristic values 
L*(L) are found from the minima of Q as a function of 
L T for a given L, and the critical coupling K c is found 
from the crossing points of these minima as a function 
of K. The correlation length exponent v is determined 
through finite-size scaling of the related quantity 



(d(m 2 )/dKY 
d(m 4 )/dK 



(3) 



where the derivatives are calculated by d(m n ) / dK = 
(E x )(m n )-(E x m n ), with E x = - cos{6 x>T -6 x+liT ). 
To extract critical exponents f3 and 7, we use the usual 
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FSS forms for the magnetization 

(|m|) oc L~^ v 
and the magnetic susceptibility 

X = LL T (m 2 ) oc L~« /h 



(4) 



(5) 



respectively. The anomalous dimension rj is then found 
from the scaling relation r; = 2 — 7/1/. We have also 
checked that the value of 77 obtained from the susceptibil- 
ity data is in correspondence with that obtained (through 
z + 77) from the critical two-point correlation function of 
the order parameter field, G(L/2) oc j^- d ~ z -n m G f 
the above observables are evaluated at L T = L* T , and we 
are careful to only use system sizes L T relatively close to 
L* in the interpolation. Using a polynomial fit of as low 
order as 3 works very well in most cases, although more 
care must be taken when extracting v. 

The error estimates we report are jackknife estimates 
of statistical errors only, but include contributions from 
the uncertainty in L* and the critical coupling K c . The 
value of K c is in general extrapolated from the scaling 
form K*(L) = K c + cL~ u ' for the crossing points K*(L) 
of Q(K,L T = L*) for adjacent system sizes L. For re- 
gions of s where the crossing points coincide and the ex- 
trapolation procedure breaks down, we base the estimate 
of K c on the largest values of L. When extracting critical 
exponents, we make sure to use system sizes large enough 
for the above mentioned FSS forms to be valid. Possible 
corrections to scaling are discussed below. We note in 
particular that we initially assume a second-order phase 
transition for all values of s we use, so critical exponents 
obtained in the case of a BKT transition should be re- 
garded as effective exponents only. The special case of 
s « 2 for the XY model is therefore re-examined sep- 
arately in Sec. 
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are discussed for a special case in Sec. 
The Monte Carlo simulations are performed using a 
Wolff cluster algorithrrPH f or long-range interactions.^ 
The results are obtained using an implementation of 
the Mersenne Twister^ random number generator, 
but other random number generators produce consistent 
results. Ferrenberg-Swendsen reweighting techniques 31 
were applied to the data. For the simulations of the XY 
case, we use a model with Z 32 symmetry to emulate the 
continuous U(l) symmetry. 



III. RESULTS AND DISCUSSION 

When extracting critical exponents for the model we 
consider in this paper, we anticipate that the only param- 
eter in Eq. relevant in determining the universality 
class is the interaction decay exponent s. (The values 
we present Monte Carlo results for are s — 1, 1.25, 1.5, 
1.625, 1.75, 1.875, and 2.) Nevertheless, we also find 
that the corrections to scaling are strongly affected by 
the strength of the dissipation term as quantified by a, 



for finite systems. In order to minimize finite-size effects 
and ensure a relatively fast onset of asymptotic values of 
the exponents, a specific value of a could in principle be 
tailored to each value of s.^ Instead of adjusting a for 
each individual value of the decay exponent s, we have di- 
vided the span of s values into two regions where we have 
applied different sets of coupling constants. For s > 1.625 
where we expect the dissipation term to be weakly rel- 
evant (in the sense of a small correction-to-scaling ex- 
ponent) or even irrelevant, we set the coupling values 
according to a = 0.1 and K T = — In (tanh 4) ss 0.7719. 
For s < 1.5 we find that it is more appropriate to choose 
a larger value of a while reducing K T in order to observe 
a rapid finite-size crossover to the asymptotic exponents. 
In this region we use a = 0.5 and K T = 0.4. For the in- 
termediate value s = 1.625, we use a = 0.3 and K T = 0.4. 
We can easily confirm for the smallest values of s that 
the universality class does not depend on the value of a, 
but corrections to scaling makes this harder for larger s, 
as discussed in Sec. IIIICI 



A. Results for the critical exponents 
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FIG. 1: Dynamical critical exponent z as a function of s for 
the Z4 and the XY model. The naive scaling estimate zo (the 
solid curve) does not coincide with the calculated z for values 
of s other than the integer-valued end points of our span of s 
values. 

In Fig. [T]we present the dynamical critical exponent z 
as a function of s. A notable feature of the results is the 
similarity between the two order parameter symmetries. 
To the accuracy of our simulations, there is essentially 
no difference between the continuous U(l) symmetry and 
the discrete Z4 symmetry. Also, the calculated z values 
do not conform to the scaling estimate zq — 2/s, but 
instead fall off faster for increasing s than expected from 
naive scaling. 

We present the evolution of 77 as a function of s in Fig. 
[2] Again, we find coinciding values for the two order pa- 
rameter symmetries. For both Z4 and U(l), rj increases 
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FIG. 2: Anomalous scaling dimension r\ as a function of s 
for the Z4 and the XY model. 770 indicates the discrepancy 
between the naive scaling estimate zo = 2/s and the actually 
calculated value of the dynamical critical exponent z (based 
on the mean for the Z4 and XY model). 



steadily with decreasing effective dimension of the sys- 
tem. Also shown in Fig. [2] is the quantity rjo = 2/s — z, 
which quantifies the difference between the naive scal- 
ing estimate zq and the calculated z. For s < 1.25 the 
evolution of 770 closely follows the calculated values of 77, 
making the scaling relation z = 2/s — rj& fair approxima- 
tion also for s > 1. For larger values of s, however, this 
scaling relation has clearly broken down, as the values of 
z again approach the naive estimate as s — > 2. In this 
limit, 77 approaches the value 77 — 1/4, which is expected 
for both the 2D Ising model at the critical point, as well 
as for the 2D XY model at the critical end point. 
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FIG. 3: Critical exponent ratio /3/u as a function of s. This 
ratio appears to be independent of order parameter symmetry 
and is also well defined in the limit of large s. The dashed 
line represents /3/u = 1/8. 

Next, we turn to the remaining critical exponents. Fig- 
ure [3]shows the results for the ratio /3/u as obtained from 
the magnetization. We do not show the ratio 7/17 al- 
though its behavior is easily inferred from Fig. [2]and the 



relation 7/1/ = 2 — 77. Again, the FSS exponent seems 
to take essentially the same values for the XY model 
as for the Z4 model. This is also the case for s — > 2, 
where we expect the dissipation to be irrelevant so that 
the effective dimensionality is reduced to d c g = 2. For 
the XY model, the U (1) symmetry of the variables then 
cannot be spontaneously broken, and the strong-coupling 
phase of the model features only quasi-long-range order 
(QLRO). Nonetheless, the system develops a finite mag- 
netization to as a finite-size effect, with a well-defined 
FSS exponent. The value (3/v » 0.125 = 1/8 of this 
exponent when s = 2 (as well as the corresponding sus- 
ceptibility ratio 7/z/ ss 7/4) is also found for the classical 
2D XY model and is, incidentally, the same as the cor- 
responding ratio in the 2D Ising model. We discuss this 
issue in more detail in Sec. IIII El 
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FIG. 4: Correlation length exponent v as a function of s. As 
the dissipation term becomes more short ranged with larger 
s, the exponent for the Z4 model approaches the 2D Ising 
value v = 1. In the XY case, v is expected to diverge in the 
limit doff — ► 2. The results (obtained for finite L) presented 
for the largest values of s should therefore be regarded only 
as effective exponents. 

The correlation length exponent v is shown in Fig. ^ 
while the critical exponent /3 is shown in Fig. [5] We 
do not show the results for the exponent 7 here, but its 
behavior is qualitatively very similar to that of the expo- 
nent v. In the Z4 case, both exponents start out close to 
the 3D Ising limit for s — 1, and approach the 2D Ising 
limit indicated by the dashed line when s — > 2. Con- 
sider now the XY case. When s = 1, these exponents 
take on values close to those of the 3D XY model. How- 
ever, the exponents /3, v and 7 are not well defined when 
d c ff = d + z = 2, as their values are formally infinite at 
a transition separating a disordered phase and a QLRO 
phase. This is the case when s = 2. Our FSS analysis for 
these exponents, which presupposes a second order phase 
transition, is strictly speaking not applicable to the BKT 
transition. The resulting (effective) exponents /3, v and 7 
appear to diverge as L — > 00 close to s — 2. Note that al- 
though £ is exponentially divergent at a BKT transition, 
we may still define z through the relation £ T oc £ z . 
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In practice, the logarithmic scaling analysis^ is usually 
best suited for small to moderate system sizes. Conse- 
quently, instead of scaling the helicity modulus directly, 
we resort to calculating k via other observables, and find 
that 2(3 /v — rj — within statistical uncertainty for s = 2. 
Moreover, 2/3 /v — rj is very close to zero for all s > 1.75. 
For even smaller values of s (where direct scaling of T 
is more reliable), we have confirmed that the scaling law 
k — 2(3 jv — 77 is valid also in the presence of dissipation. 
This scaling form is also equivalent to k = d c g — 2. Thus, 
the helicity modulus vanishes continuously as K — > K([ , 
provided d e ff = d + z > 2. The above equivalence as- 
sumes that hyperscaling is valid, and we have confirmed 
this validity for all values of s. 



FIG. 5: Critical exponent /3 as a function of s. For the Z4 
model, /3 evolves smoothly from the 3D Ising to the 2D Ising 
limit as s is increased from 1 to 2. The XY result, on the 
other hand, starts out near the 3D XY value for s = 1 and 
features a nonmonotic evolution of j3 with s, with a divergent 
(S in the limit of large s; cf. Fig. [4] 

Another observation in the U(l) case is that while 
the combination (3/v is monotonically decreasing with 
increasing s, f3 itself is exhibiting a non-monotonic evo- 
lution as a function of s. The value of j3 is at first decreas- 
ing as the increasing value of s drives the system away 
from the 3D behavior, just as for the Z4 case. However, 
as mentioned above, /3 is divergent in the 2D XY limit, 
and the reduction of (3 is therefore reversed at an inter- 
mediate value of s 



B. Berezinskii-Kosterlitz-Thouless phase transition 
and the helicity modulus 

The helicity modulus 

T = zrfe cos( ^" 9iHr) ) (6) 

T \x=1t=1 / 

is expected to scale as T oc L K at a critical point 
where a U(l) symmetry is spontaneously broken, with 
k = 2(3 /v — r\. For a 2D XY model, however, the expo- 
nent k is exactly zero, reflecting the fact that at the BKT 
phase transition, the helicity modulus jumps to a finite 
value with logarithmic finite-size corrections. By direct 
comparison of the calculated values of T for s = 2 and 
the scaling form expected for a BKT transition,^ un- 
ambiguous conclusions regarding the universality class of 
the phase transition at s = 2 could not be drawn. The 
presence of the presumably irrelevant dissipation term 
is still effective in driving the system away from BKT- 
type criticality at all but the very largest system sizes. 



C. Boundary between long-range and short-range 
critical behavior 

From Figs. [T] to [5] it is evident that all critical ex- 
ponents are very close to their short-range values for 
s > 1.75. The naive scaling estimate places the boundary 
at which the dissipation term becomes irrelevant at s = 2. 
For classical models with (isotropic) long-range interac- 
tions d ecayin g with distance r as l/r 1+s , it has long been 
debatecP^EH whether the models feature the exponents of 
the corresponding short-range model already when s ex- 
ceeds a value s* = 2 — r)sn- Here, denotes the anoma- 
lous dimension of the short-range model. Using large- 
scale Monte Carlo simulations, it has been showrPH that 
for the long-range 2D Ising model, the anomalous dimen- 
sion follows the conjectured exacipS relation 77 = 2 — s for 
s < 1.75, but that r) = 0.25 = 77SR for s > s* = 1.75. Al- 
though the long-range interaction of the dissipativc quan- 
tum models we consider is highly anisotropic, in contrast 
to the isotropic classical long-ranged models, it is plau- 
sible that also in these models the threshold value of s 
beyond which dissipation is irrelevant is reduced from 
s = 2 to some lower value. 

In order to establish this boundary more accurately 
also for the present case of anisotropic interactions, we 
have performed a more careful analysis of the case s = 
1.875 for the Z4 model. Including corrections to scal- 
ing, using the ansatz L* = aL z (l + bL^^), we find 
z = 1.002(11). Hence, the decay exponent s — 1.875 
may serve as an upper bound for the boundary value s* 
necessary to render the dissipation term effectively short- 
ranged. This, in turn, would render the system effec- 
tively Lorentz invariant with z = 1. To get the statistics 
required to include corrections to scaling in a stable man- 
ner, we included three different values of the dissipation 
strength a in a joint fit. This also provides an a posteriori 
justification of the choice of lower values of a for higher 
values of sP^ Probably due to logarithmic corrections 
expected at the presumed boundary value s* = 1.75,^ 
we are not able to acquire the same level of accuracy 
for this value of s. Therefore we cannot rule out that 
the dissipation term is rendered effectively short ranged 
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at some other value s* £ (1.75,1.875). An exceedingly 
slow crossover to asymptotic critical exponents for values 
s ~ s* can conceivably be understood from the competi- 
tion between the fixed point corresponding to short-range 
(Lorentz-invariant) criticality and the fixed point corre- 
sponding to long-range (dissipation-dominated) critical- 
ity. 

We close this section with a remark on the evolution 
of the anomalous dimension. In the quantum dissipa- 
tive model we have studied, the anomalous dimension 
increases for increasing s. This is a consequence of the 
effective dimensionality d c g = d + z decreasing with in- 
creasing s. Lowering the dimensionality from the upper 
critical dimension, where 77 = 0, tends to increase r\. This 
is quite different from the situation encountered in clas- 
sical models with isotropic long-range interactions. Clas- 
sical models with short-range interactions and an action 
of the form S ~ q 2 <fi q (f>- q (where <j) q is an appropriate 
order-parameter field) have propagators G(q) ~ 1 / 1 <y | 2 77 . 
The corresponding long-range models with an action of 
the form S ~ \q\ s <j) q 4>- q have propagators G(q) ~ 1/M S , 
when s < 2 — rygR,. One may now, as is customarily done 
in the literature on long-range classical isotropic models, 
define an effective anomalous scaling dimension for such 
systems by comparing with the corresponding expression 
for the short-range case, finding r\ = 2 — s, which de- 
creases with increasing s. This relation is best viewed 
as a result of somewhat artificially imposing the stan- 
dard scaling form of a propagator for short-range systems 
(l/|g| 2- '') on the form of the propagator for systems with 
long-range interactions, 1/|<7| S . 



D. Scaling relation between z and 77 

In Fig. [2j we demonstrated how the scaling relation 
z — zq — r] cannot be valid except close to the Ohmic 
limit s = 1 and that 77 ss 1/4 for all s > 1.75. Moreover, 
from our numerics we think it is likely that z(s) = 1 for 
s > 1.75. A scaling relation between z and 77 which would 
fit well with these observations is 
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The scaling relation z — (2 — rj)/s has been suggested 
previously in Ref. |T7] in the context of a damped non- 
linear a model. In Fig. [6j we show the same data for 
the dynamical exponent z as in Fig. [I] but compared 
with the ansatz ([7]) instead of the naive scaling estimate 
zq = 2/ s. Although there are probably still some correc- 
tions to finite-size scaling, Eq. ^ seems to fit the data 
far better than the alternatives. 

We next provide a heuristic argument for why the scal- 
ing relation ^ may be reasonable for s < 2 — 773^. Build- 
ing on the arguments in Sec. Ill C for classical isotropic 
long-range-interacting systems, we take as a starting 
point that a dissipative quantum model with action of 
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FIG. 6: The dynamical critical exponent 2 as in Fig. [T] but 
compared with the scaling estimate z' = (2 — rj)/s. The values 
of r\ used in the scaling estimates are the same as reported in 
Fig. [2] with the left panel corresponding to the Z4 model and 
the right panel to the XY model. 



anisotropic long-range-interacting system. Introducing 
the suitably chosen frequency coordinate uj s = q 2 + \co\ s , 
the propagator takes the isotropic form G(oj) ~ l/|w| s . 
Recall that in the quantum case, the anomalous scaling 
dimension 77 is defined from the spatial correlation func- 
tion G(x) ~ i/j.d+z-2+r] = \/ x Bx . To find 77, we Fourier 
transform the propagator to obtain the imaginary-time 
correlation function. In terms of the frequency coordi- 
nate u>, the system effectively has d' eS = 1 + d/z di- 
mensions. Therefore, the correlation function decay ex- 
ponent in terms of the isotropic space-time coordinate 
f = (t 2 + x 27 -) 1 / 2 would be 8f = (z + d)/z — s. Compar- 
ing with the imaginary-time decay exponent 8 T = 9 x /z — 
(d + z — 2 + r])/z, we find sz = 2 — 77, which is equivalent 
toEq. 0. 

Finally, we point out that our results for the scal- 
ing relation ^ for dissipative models should also have 
relevance for nondissipative q uantu m spin chains with 
long-range spatial interactions! 38 * 3 ^ One arrives at ex- 
actly such a model by simply interchanging the x and 
r coordinates of the action we have considered. The 
dynamical critical exponent of this model is given by 
z' = l/z = s/(2 — 77), with the quantity r/(s) evolving 
as shown in Fig. [2] This quantity will, however, not 
be identical to the anomalous scaling dimension of the 
model, 77', which is given by the classical result 77' = 2 — s. 
Hence, there exists no independent scaling relation be- 
tween the dynamical critical exponent z' and the anoma- 
lous dimension 77' for a nondissipative quantum spin chain 
with long-range interactions. 



the form S ~ (q 



s )4'q,ui4 l -q,-u) can be viewed as an 



E. Dependence on symmetry 

For s — 1, we have z w 2, d e ff ~ 3, while for s = 2, 
we have z = 1, d c g = 2. For these two cases, it is known 
either analytically or numerically that the exponent 77 is 
very similar for the Ising and XY modelsPS There ap- 
pears to be no particular deep reason for this. For in- 
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stance, the well-known value 77 = 1/4 comes about for 
completely different reasons in the 2D Ising and 2D XY 
models, and their similarity thus appears to be acciden- 
tal. Using the scaling relations 2/3/f — rj = d e g — 2 [as- 
suming Eq. ^fj] and j/u = 2 — rj, it follows that the 
similarities in j3/v and 7/1/ for the Ising and XY models 
are as coincidental as they are for rj, both in 2D and 3D. 
It appears that these coincidences persist in all dimen- 
sions between 2 and 3. There is good reason to expect 
that the same also holds in the sub-Ohmic regime s < 1. 
Such values of s increase the effective dimensionality be- 
yond 3, eventually driving all exponents to their universal 
mean-field values at the upper critical dimension. 

We next comment on other values of q, and how our 
results apply to those cases. The Ising and XY models 
represent extreme cases of Z q clock models, with q = 2, 
q = 00, respectively. The partition function for the q = 4 
case is simply the square of the case q — 2, and hence 
they give identical results. For larger q > 4, anisotropy 
is irrelevant and we thus expect the results of U(l) to 
emerge. We therefore conjecture that the results of this 
paper for z, rj, [3/v, and 7/1/, are valid for all Z q clock 
models. The only possible exception is the case q = 3, 
also equivalent to the three-state Potts model, where the 
anisotropy with respect to a U(l)-symmetric model is 
known to be relevant. Although we have not checked 
this, it may still be possible that Eq. Q holds also for a 
dissipative Z3 clock model, at least for s > lp2 

An alternative perspective on this, supporting the no- 
tion that the scaling relation z — (2 — r/)/s is valid for all 
q, may be provided by the following qualitative argument. 
The variation of z with the parameter s determining the 
range of the dissipation, expresses a variation in the ef- 
fective space-time dimensionality of the system. This is 
determined by the interaction of the spins at a given site 
in the imaginary-time direction. Due to the long-range 
character of this interaction, each spin interacts with a 
large number of fluctuating copies of itself along a chain 
in the imaginary-time direction. Due to the summation 
over many spins at different Trotter slices, the discrete 



nature of the spins in a Z q clock model is washed out, 
even in the case q — 2. Therefore, the manner in which 
the dissipation affects the effective dimensionality does 
not depend on whether the spins at each space-time lat- 
tice point take on discrete or continuous values. 



IV. SUMMARY 

We have performed Monte Carlo simulations on a gen- 
eralized spin-boson model in one spatial dimension fea- 
turing non-Ohmic site dissipation and two variants of 
order parameter symmetry, namely Ising- like and U(l). 
By tuning the imaginary-time decay exponent of the dis- 
sipative interaction, s G [1,2], we are able to contin- 
uously vary the effective dimensionality of the system. 
Apparently, the order parameter symmetry has very lit- 
tle bearing on the evolution of the effective dimensional- 
ity, d e f[ = d + z, as a function of the decay parameter s. 
While naive scaling estimates a crossover from criticality 
dominated by the dissipation term to an irrelevant dissi- 
pation term at s = 2, we measure exponents in relatively 
good correspondence with the underlying, short-range in- 
teracting model at a somewhat lower value s « 1.75. Our 
results also suggest that for 1 < s < 2, the exponents z 
and rj to a good approximation obey the scaling relation 
z = max {(2 — r/)/s, 1}. 
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